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ABSTRACT 

Graham et al. have detected a 5.2 yr periodic optical variability of the quasar PG 1302-102 at 
redshift 2 : = 0.3, which they interpret as the redshifted orbital period (1 -|- 2:)tbin of a putative 
supermassive black hole binary (SMBHB). Here, we consider the implications of a 3 — 8 
times shorter orbital period, suggested by hydrodynamical simulations of circumbinary discs 
(CBDs) with nearly equal-mass SMBHBs (q = M 2 /Mi > 0.3). With the corresponding 2 —4 
times tighter binary separation, PG 1302 would be undergoing gravitational wave dominated 
inspiral, and serve as a proof that the BHs can be fuelled and produce bright emission even 
in this late stage of the merger. The expected fraction of binaries with the shorter fbin, among 
bright quasars, would be reduced by one to two orders of magnitude, compared to the 5.2 yr 
period, in better agreement with the rarity of candidates reported by Graham et al. Finally, 
shorter periods would imply higher binary speeds, possibly imprinting periodicity on the light 
curves from relativistic beaming, as well as measurable relativistic effects on the Fe K a line. 
The CBD model predicts additional periodic variability on time-scales of tbin and ^ 0.5tbin, 
as well as periodic variation of broad line widths and offsets relative to the narrow lines, which 
are consistent with the observations. Future observations will be able to test these predictions 
and hence the binary-i-CBD hypothesis for PG 1302. 
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1 INTRODUCTION 

[Graham et al.| ( |2015] hereafter G15) recently reported strong optical 
variability of the quasar PG 1302-102, with an observed period of 
^obs = 5.2zb0.2yr. G15 attribute the variability to the orbital mo¬ 
tion of a super-massive black hole binary (SMBHB). Broad emis¬ 
sion lines in the spectrum of PG 1302 imply a binary mass in the 
range M = Assuming that the binary’s orbital pe¬ 

riod tbin equals the rest-frame optical variability period topt, G15 
derive a fiducial binary separation a ^ (0.0084 zb 0.0003)pc 
(276 zb 9)Rs for M = lO^'^M©, where Rs = ^GMjc^ is the 
Schwarzschild radius. 

Hydrodynamical simulations of a binary BH embedded in a 
gaseous accretion disc predict that, depending on the binary mass 
ratio and the physical parameters of the disc, the strongest periodic¬ 
ity in the accretion rate on to the BHs may correspond to the motion 
of gas farther out in the disc, at a few times the binary separation, 
producing optical variability at several (~3-8) times the binary or¬ 
bital period. In this article, we discuss this expectation, and show 
that a reduced binary period, in the case of PG 1302, would have 
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several important implications. Follow up spectroscopy and photo¬ 
metric monitoring can determine the true binary period. 

The rich variability structure of the mass accretion rates seen 
in simulations can be roughly divided into four distinct categories, 
based on the binary mass ratio q = M 2 /Mi . For q <0.05, the disc 
is steady and the BH accretion rate displays no strong variability 
( jD’Orazio et al. [2013] [Farris et al.|2014| D’Orazio et al. in prepa¬ 
ration). For 0.05 < ^ < 0.3, the accretion rate varies periodically 
on the time-scale tbin, with additional periodicity at 0.5tbin. Bi¬ 
naries with 0.3 < q < 0.8 clear a lopsided central cavity in the 
disc, causing variability on three time-scales. The dominant period, 
(3 — 8)tbin is that of an over dense lump, orbiting at the ridge of 


the cavity, with additional periodicities at tbin and 0.5tbin 

([Mac-[ 

Fadyen & Milosavljevic|2008|[Shi et al.|2012| Noble et al. 

2012| 

Roedig et al.|2012| [D’Orazio et al.|2013| [Farris et al.|2014 

1 . The 


dominant period depends on the size of the cavity, and thus on disc 
parameters, such as temperature and viscosity. Finally, equal-mass 
(g = 1) binaries display variability at the longer lump period and 
at ~ 0.5tbin' 

Here, we consider the identification of the observed variability 
of PG 1302 with the long, cavity-wall period, and introduce the 
parameter x = ^opt/^bin, denoting the ratio, 3 < x ^ 8, of the 
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Figure 1. Results of 2D hydrodynamical simulations of a binary BH sur¬ 
rounded by a circumbinary accretion disc. The BHs clear out a central cav¬ 
ity and form their own minidiscs. Left-hand panels: snapshots of the (log¬ 
arithmic) surface-density of the gas discs, after reaching quasi-steady state, 
with mass ratios of q = 0.5 (top) and q = 0.1 (bottom). Right-hand pan¬ 
els: corresponding LSPs of the total accretion rate on to the secondary + 
primary BHs. The discs are locally isothermal with a Mach number of 10 
and an alpha viscosity prescription (a = 0.1). 


observed rest-frame period and the true binary period. The binary 
separation is then 

/Y\-2/3/ M 

.»(SI4±3)ife(|) , <„ 

or (0.0029 zb 0.0001) pc for the fiducial choices of x and M. 

In the rest of this article, we first (^ explore several impli¬ 
cations of a reduction in the binary’s orbital period, including the 
nature of PG 1302’s orbital decay and its ability to produce elec¬ 
tromagnetic (EM) radiation ( g2.1| ), the expected binary fraction of 
quasars ( ^2.2| ), and the detectability of gravitational waves (GWs) 
from PG 1302 by pulsar timing arrays (PTAs). We then (^ pro¬ 
pose possible observational tests of the underlying binary BH -i 
circumbinary disc (CBD) model, including variations of broad line 
widths and centroids correlating with the optical variability ( ^3.1| ), 
additional periodic variability at the true tbin caused by relativis¬ 
tic beaming ( ^3.2| ), signatures in the broad Fe K a lines ( ^3.2| ), and 
the existence of distinct secondary peaks in the periodogram ( g3.3| ). 
We briefly summarize our main conclusions in 0 


wider range of simulations is needed in the future, to address possi¬ 
bilities such as eccentric ( [Roedig et al.|2012| ), tilted ( [Hayasaki et al.| 
|2015| ), or retrograde { [Nixon et al.|2011| ) binary orbits. 

The results are illustrated in Fig[2 The left-hand panels show 
snapshots of the surface density, and the right-hand panels show 
Lomb Scargle periodograms (LSPs) of the total accretion rate mea¬ 
sured in the two BH minidiscs over 200 binary orbits. The top pan¬ 
els, for g = 0.5, show an over dense lump orbiting at the rim of 
the central cavity, resulting in strong periodicity at the orbital time 
~ Gtbin of the lump. The periodogram shows weaker peaks at tbin 
and at ~ O.Gtbin- This three-time-scale behaviour, with the longest 
time-scale dominating, is observed for 0.3 <q < 0.8. The location 
of the highest frequency peak is closer to 0.5tbin near the low end 
of this range (q ~ 0.3), and also has a weak dependence on disc 
temperature and viscosity which must be quantified in future work. 
The bottom panels, for = 0.1, show no orbiting lump and exhibit 
accretion rate periodicity only at tbin and 0.5tbin. This behaviour 
is found in the range 0.05 <^<0-3. 

[Farris et ^ ( |2014j l have shown that for unequal-mass binaries, 
accretion occurs preferentially on to the secondary BH, with the 
ratio of accretion rates as skewed as M 2 /Mi 10 — 20 in the 

range 0.02 < q <0.1. Over long time-scales, this would drive 
the binary to more equal masses, suggesting that mass ratios of 
0.3 < q < 0.8 may be common. Near-equal mass binaries are 
also preferred in cosmological models of the population of merging 
SMBHs ( [Volonteri et al.[2003| ). 

Although there are large uncertainties in how accretion rate 
fiuctuations turn into luminosity variations, we do expect the sim¬ 
ulated accretion rate variations to lead to optical luminosity varia¬ 
tions for PG 1302. The luminosity will follow local accretion rate 
fiuctuations when the longer of the thermal or photon diffusion 
time-scale is much shorter than the accretion rate fiuctuation time- 
scale (~ tbin). This is indeed the case where accretion modula¬ 
tions occur in our simulations, at the minidisc edges. Furthermore, 
optical emission is generated at the minidiscs edges. We compute 
thin disc spectra for the CBD and minidiscs. For the preferred mass 
range of PG 1302, near unity mass ratios, and the expected range of 
X, the dominant optical component of the spectrum is generated by 
the low-energy tail of the blackbody emission from the outer edges 
of the minidiscs, as well as (steady) emission from the inner re¬ 
gions of the CBD. Disc-fbinary simulations by |Farris et al.| { [2015b] ), 
which self-consistently compute the local effective disc tempera¬ 
ture not assuming a steady state, find results in agreement with 
our analytic reasoning: luminosity variations track the accretion- 
rate fluctuations, except at low frequencies where the quasi-steady 
CBD dominates. 

The above lines of evidence motivate us to examine the possi¬ 
bility that the apparent 5.2 yr period in PG 1302 is the (redshifted) 
lump period, and assess the implications. 


2 IMPLICATIONS OF A SHORTER ORBITAL PERIOD 

In order to demonstrate the possibility of a short orbital period for 
the PG 1302 binary, we have performed hydrodynamical simula¬ 
tions, following the set-up in our earlier work ( [Farris et al.|2014| ). 
The hydrodynamical equations are evolved for > GOO binary orbits, 
using the two-dimensional code DISCO { [Duffell & MacFadyen| 
[2011| ), with the BHs moving on fixed circular orbits, surrounded 
by an isothermal (Mach number =10) disc, obeying an a-viscosity 
prescription (a = 0.1). The fluid motion around individual BHs is 
well resolved (with a log grid of 384 radial cells extending to 8a 
and a maximum of 512 azimuthal cells). The two runs discussed 
below differ only in their BH mass ratio (q = 0.1 and q = 0.5). A 


2.1 Binary-Disc Decoupling 

A shorter orbital period would place the binary at a later stage 
of its orbital decay. A critical point during the orbital decay is 
the decoupling of the binary from the CBD, and it is important 
to know whether the binary is past this point. Here, we consider 
the decoupling radius for which the GW decay time-sc ale becomes 
shorter than the decay time-scale due to gaseous torques (so-called 
secondary-dominated Type II migration; |Syer & Clarke|19^ , out¬ 
pacing the CBD. We use simple ID models of the binary -1 disc 
system ( [Haiman et al.[2009[ hereafter HKM09) to calculate the sep¬ 
aration row at which decoupling occurs for a circular binary with 
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Figure 2. The ratio of binary separation a to the decoupling radius row? 
for three different values of the ratio between the rest frame optical period 
and the true binary period, x = fopt /f bin • The shaded region marks the 
binary mass range inferred from the widths of broad lines measured by G15. 
For X > 3, the PG 1302 binary is past decoupling, for any choice of mass 
M. 


mass ratio q = 0.3. We assume an ct-viscosity ia = ctPgas(p^)~^, 
with gas pressure Pgas, density p, and disc angular velocity Q. All 
other disc parameters are assumed to have the fiducial values given 
in HKM09. In Fig.|^ we plot the ratio a/row as a function of the 
total binary mass M, with the binary separation from equation Q. 
for the range of masses in G15, and for three values of x covering 
the range suggested by the hydrodynamical simulations. 

Interpreting the observed variability in PG 1302 with tbin, as 
may be justified for q <0.3, it is unclear whether or not the binary 
has entered the GW dominated regime and decoupled from the disc. 
The binary would still be coupled to the disc for M < 10®'^M© 
(the majority of the range inferred from the broad lines by G15), 
but GW-driven and decoupled if M > lO^'^M©. However, for the 
shorter binary periods 3 < x ^8, justified for 0.3 < q <0.8, 
we find that the inferred smaller binary separation would place the 
binary well past decoupling. For g > 0.3 and a < 0.3, the binary 
in PG 1302 is plunged even deeper into the GW-dominated regime. 

Because the binary outpaces the disc, it has been argued that 


the post-decoupling BHs may be “sta rved” and thus dim (|Milosavl-| 


[jevic & Phinney|2005|[Shapiro|2010 

[Tanaka & Menou 

2010[l. Re- 

cent simulations |Noble et al.|2012| 

Parris et al.|2015a 

1 show that 


high levels of accretion can persist well past the decoupling phase, 
delivering gas to the binary efficiently until much closer to coa¬ 
lescence. These simulations also exhibit the lopsided cavity which 
generates the x^bin variability considered here. Identification of the 
variability in PG 1302 with the cavity wall lump period would con¬ 
stitute the (to our knowledge, first-ever) detection of an SMBHB 
which is undergoing GW dominated inspiral, yet producing bright 
emission, near the Eddington limit[^ 

2.2 Binary Fraction among Quasars 

A shorter binary orbit also reduces the expected number of de¬ 
tectable SMBHBs in quasars, because gas driven binaries are ex¬ 
pected to spend less time at smaller separations. A simple estimate 
of the fraction of quasars that would harbour binaries with an or¬ 
bital period tbin can be obtained from the residence time Ues = 
a/d at each separation a, and the lifetime of bright (Lg/Ledd > 


0.3) quasars, tq ~ 10® yr (|Martini|2004|l. Assuming that a fraction 
/bin of all quasars are triggered by coalescing SMBHBs (e.g |Hop-| 
[kins et al.|2007| and references therein), it follows that among bright 
quasars, the fraction with orbital period tbin is /var — /bin/duty 5 
where /duty = Aes(tbin)/tQ is the fraction of the bright quasar 
phase that a typical binary quasar spends at the orbital period tbin- 

We use the binary-i-disc models of HKM09 to predict the res¬ 
idence times of binaries. Prior to decoupling, tres is determined by 
the binary’s interaction with the gas disc. For the masses and sepa¬ 
rations relevant for PG 1302, the disc would be radiation pressure 
dominated, yielding a relatively shallow power-law dependence 
tres oc with 0.5 < ^ < 1.5. These scalings depend on the 
poorly understood physical model of the disc and its coupling to 
the binary. Past decoupling, the residence time is precisely known, 
since it is determined by the strength of GWs. The dependence is 
much steeper, tres oc t^.^^ oc ■ For reference, a binary with 

M = lO^'^M© and tbin = 4yr would be in the disc-driven stage 
and would have tres ~ 10® yr, yielding a large /duty ~ 10“^. 

The expected /var can be compared with the number of peri- 
odic candidates uncovered in CRT S ( [Drake et al.|2009[[Djorgovs]a| 
et al. 12011 Mahabal et al. |20 if] ) p] There are 114,000 quasars 


in the CRTS sample with luminosity higher than PG 1302, 6 

of which are SMBHB candidates with period tobs < 5yr (Gra¬ 
ham, private communication), amounting to an observed fraction 
of f?ar = 5 X 10“®. Fig. [^illustrates combinations of M, x, and 
/bin, for which we expect 1 (light grey regions) or 10 (dark grey) 
candidates in the CRTS quasar sample with periods ^ 5.2 yr. Each 
shaded region is bounded by the assumed fraction of quasars re¬ 
lated to SMBHBs at all, /bin = 0.01 (left) to /bin = 1 (right). 

If the observed period of PG 1302 is assumed to be the binary 
orbital period (x = 1), then the rarity of the binary candidates in 
CRTS require /bin < 0.14 (< 0.19) at g = 0.3 (q = 1.0), even at 
the most extreme mass of M = 10® ^M©. Taking the G15 fiducial 
mass of M = 10®'^M©, these fractions must be as low as /bin < 
0.006. These low values would be surprising, as a large fraction of 
quasars are commonly believed to be triggered by mergers. This 
association is based on various pieces of observational evidence, as 
well as on the success of merger-based quasar population models to 
reproduce many properties of the observed quasar population (e.g. 
[Kauffmann & Haehnelt|2000] |. If, instead, the observed period of 
PG 1302 is due to the 3 — 8 times longer lump-periodicity, then 
the SMBHB fraction and the inferred binary mass of PG 1302-12 
come into wider agreement with the expectation that /bin = 0(1). 
e.g. allowing /bin ~ 0.3 with q = 0.3 and M = 10® ®M©. 


2.3 Detectability of GWs 

A reduced orbital period increases the frequency and amplitude 
of GWs emitted by a binary, and it is interesting to ask whether 
PG 1302 may be detectable by present or future PTAs. The 
GW frequency, /gw = 61 (x/5) nHz, places the bi¬ 

nary in the range of PTA sensitivity (e.g. [Hobbs et al.1[M0l l. 
We calculate an SNR for the PG 1302 binary choosing opti¬ 
mistic binary parameters M = 10®‘^M©, q = 0.5, and tbin = 
1447x~^ d. The GW induced rms timing residual (for simplic¬ 
ity, adopting the sky and polarization averaged values) is Stow = 


^ Note that in the precessing binary model for OJ287 jLehto & Valtonen[ 
[1996[[Wltonen et al.[2008) , the orbital period is 12.2 yr, the primary is very 
massive (~ 1.8 x 10^^ M©), but the secondary is light (~ 1.4 x 10® M©). 
The latter reduces the efficiency of GWs, but increases the impact of a gas 
disc; as a result, the OJ287 binary is gas-driven, well before decoupling. 


^ Since the expected /var(Af, fbin) declines steeply with increasing M 
and decreasing fbin, it is a good proxy for the fraction of quasars with 
period fbin or less, and BH mass M or higher (or equivalently luminosity 
L or higher, further assuming a monotonic relation between L and M). 
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Figure 3. Combinations of total binary mass M and x = fopt /f bin for 
which the predicted binary fraction of quasars matches 10 CRTS candi¬ 
dates with luminosity above that of PG 1302 (dark grey, bounded by solid 
curves), and for which it consists only of PG 1302 (light grey, bounded by 
red dashed curves). Vertical lines delineate the range of masses preferred by 
broad line widths. Each shaded region is bounded by the fraction /bin of 
quasars which are triggered by a binary. In each case, the lines correspond 
to /bin = 0.01, 0.1,1 (left to right). 


V^Vi5V(27t/Gw)^/7G^ W 2.6 (x/5)^/® ns for a one 
year observation time Tyr; for /gw >10 nHz, the timing residual 
noise is nearly constant in frequency. Using noise curves of cur¬ 
rently operating PTAs, the SNR for GW detection of the PG 1 302 


binary is ~ 0.005 (x/5)^^® for NANOGrav (fig. 12 in 


Ar- 


zoumanian et al. [2014 1 or ~ 0.011 PPTA 


(black curve m fig. 9 oi |Zhu et al.|2014| ). The reduced binary period 
increases the SNR, but this increase is unfortunately modest. Fu¬ 
ture detectors, such as the international pulsar timing array (iPTA; 
[Manchester & IPTA|2013|| as well as inc lusion of the square kilo- 
metre array (SKA; |Dewdney et al.|20()^ in the PTA telescope net¬ 
works will improve the SNR by about an order of magnitude, but 
PG 1302 will remain a factor of ~ 10 below detection. 


3 TESTING THE BINARY BH SCENARIO FOR PG 1302 

3.1 Broad Line Variability and Asymmetry 

report a ~ (150 zb 50)km s“^ offset between 
the broad and narrow line components of PG1302’s H/3 emission 
line. This is much smaller than the secondary BH’s orbital speed, 
V2 = 14, 500 (1.5/[1 + g])(M/10®-W0)i/yx/5)^^®km s-\ or 
the width of the broad lines. Such larger offsets have been predicted 
for binary SMBHs, assuming that the broad line region (BLR) orig¬ 
inates from gas bound to one component of a binary SMBH and 
thus shares its overall orbital motion (e.g jTsalmantza et al.|2011j l. 
Here we argue that the smaller offset for PG1302 can also be at¬ 
tributed to a binary SMBH, assuming that the BLR is located far¬ 
ther out, in the circumbinary gas. Using simple toy models, we 
show that the lopsided geometry of the CBD gas (Fig.[^ could gen¬ 
erate the small observed offset but large width of PG 1302’s broad 
H/3 line. The models below are meant to be mere illustrations; a 
self-consistent description of the BLR is left to future work. 

The idea is that the large width of a line can refiect the orbital 
speed of gas in the CBD (over a range of annuli), whereas the offset 
of the line centroid is caused only by departures from axisymme- 
try and can be much smaller. (In a strictly axisymmetric BLR, the 
blue- and redshifts from gas on opposing sides of the binary would 
be the same and leave no net offset). To illustrate this, we compute 
the line offset Vo as the emission-weighted line-of-sight (l.o.s.) 
velocity. 


\1992\ 


Jackson et al. 


V Itz rdrdcf) 

V Ik rdrdrf) 


and the width F as the weighted rms l.o.s. velocity 

/ r \ ^ _ V P^{v<l>/rT (^io3 - Vof rdrd4> 

\2V21n2/ rdrdij} 

over a surface patch 7^ in a = 0.5 simulation (see Fig.[^. 

In order to select the scaling indices n and m, and to iden¬ 
tify a patch IZ corresponding to the BLR, we first consider a steady 
thin disc model for the CBD (as in and assume that the CBD 
is illuminated by a central ionizing source (i.e. the minidiscs). For 
PG1302’s parameters, the inner region of the CBD would have den¬ 
sity n ~ and would be highly opaque to ionizing ra¬ 

diation. In this case, recombinations in a volume corresponding to a 
very thin (AR R) inner annulus of the CBD would balance the 
central ionizing photon rate. In particular, the disc would absorb the 
covering fraction 2nRH/{4nR‘^) = 0.0bH/{0.1R) of the central 
luminosity. PG 1302’s bolometric luminosity is ~ 6x 10^® erg s 
assuming that > l/3rd of this is emitted in the UV, > 5% of which 
is absorbed by the CBD (i.e. H/R > 0.1), this would be sufficient 
to provide the total power ~ 10^^ erg s“^ measured in the broad 
lines ( jWang et al.|2003j ). The line emission from each patch of the 
CBD depends only upon the number of ionizing photons incident 
on the CBD in that direction, i.e. proportional to the scale height H 
of the inner wall of the CBD. For an adiabatic scale height, assum¬ 
ing vertical disc hydrostatic equilibrium, n = 1/3 and m — —1. 
To be specific, IZ is chosen by excising the binary plus minidiscs, 
and imposing a surface density range H/Ho = (0.01,0.5). While 
somewhat ad hoc, we find that this E range accurately picks out 
the streams edges inside the cavity, and the thin inner edges of the 
CBD. 

Although a standard Shakura-Sunyaev CBD is optically thick 
outside of the inner cavity, we also consider, for generality, an al¬ 
ternative scenario, where the BLR emission is produced by an op¬ 
tically thin medium, but still resembling the lopsided geometry in 
our simulations. In this case, line emission would scale with the re¬ 
combination rate, with n — 2 and m = 0 in equations.We adopt 
the region IZ to be an annulus with inner and outer radii (2a, 6a), 
chosen to encompass the inner CBD. 

We use the simulated surface density E, azimuthal velocity 
V(f) and the inferred isothermal scale height of the disc H = rc^jvf^ 
(assuming vertical hydrostatic equilibrium), to compute the volume 
density p — 12 jH and l.o.s. velocity uios = VcfyCoscj). All l.o.s. 
velocities are multiplied by an additional factor of sin i, where i is 
the CBD inclination angle, measured from face-on. We calculate Vo 
and r 10 times per orbit for 20 orbits. Fig. [^displays the variations 
of Vo and F with time for the optically thin case. The line centroid 
varies with mean and range 


^thn 

^thk 


( 2 zb 167) 
(-5 zb 285) 


sin i 

sin(14.1°) 
sin i 

sin(11.40) 


km / M ^ ^ ^ ^ 

V \94Rs) 


while the line full width at half-maximum (FWHM) fiuctuates pe¬ 
riodically on the lump’s orbital time. 


r*n = (4450 ±615) [^4^^ 
r*k = (4450 ±636) 


km / M \ 
\^108-5Mo ) 



The fiducial inclination angles above are chosen to match the ob¬ 
served H/3 FWHM. 
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Figure 4. Predicted variations of the centroid Vq (left) and FWHM F (right) 
of an emission line emanating from the inner CBD. The total accretion rate 
on to both black holes is over plotted in the right-hand panel in arbitrary 
units (orange). Dark black lines are smoothed versions of the light grey 
simulation data. 

We find that both models of a CBD BLR predict a line offset 
consistent with that observed in [Jackson et al.| ( |1992| ) and which 
require a CBD inclination angle that also predicts consistent line 
widths. Additionally, we find that L varies by up to ~ 14 per cent 
of the mean in each case. It is important to emphasize that in the 
optically thick case, these results arise from the non-axisymmetric 
velocities in the gas that trace the lopsided inner wall of the CBD, 
whereas in the optically thin case, they are driven by the lopsided 
density distribution. 

Since fiuctuations in the latter case arise from the lump’s vary¬ 
ing position along the cavity wall, they correlate with long-term 
variations in the BH accretion rate. In the right-hand panel of Fig. 
1^ we plot the accretion rate on to the BHs, together with the F 
variations of the optically thin case. The phase lag between line- 
width maximum and accretion maximum derives from the time 
between lump passage near the BHs, and the lump-enhanced ac¬ 
cretion in the minidiscs. The phase difference is therefore indepen¬ 
dent of the observer’s viewing angle. This is not true for the ampli¬ 
tude and shape of the FWHM variations, and the relative phase of 
Vb(t), which depend on viewing angle. For the optically thick case, 
the FWHM variations are similar, except that they have a ~ 20% 
higher mean, and are ~half a cycle out of phase with the accretion 
rate modulation. The phase difference arises because the optically 
thick model tracks the low density gas at the cavity and stream 
edges instead of the high density lump. The centroid variations are 
also similar in the optically thick case but have a few times higher 
mean and deviation because they track the eccentric cavity shape 
rather than a symmetric annulus. 

The line characteristics which we calculate here are dependent 
on the existence and magnitude of the orbiting cavity wall over¬ 
density (requiring 3 < x ^ 8). As long as the binary has mass ratio 
such that a cavity wall lump is generated, our BLR calculation is 
largely unchanged. In addition to mass ratio, disc viscosity and tem¬ 
perature affect the lump size and thus the magnitude of broad line 
variations. The line shape also depends on the broad line emission 
model. Hence, a full study of CBD broad lines, which examines 
more sophisticated recombination models and a range of disc pa¬ 
rameters, binary mass ratios and viewing angles, is warranted in a 
future study. Parameter dependences aside, observation of line vari¬ 
ability, matched to luminosity variability, would provide evidence 
for the CBD model and the origin of the BLR as well as identifica¬ 
tion of the CBD cavity wall period. 

3.2 Relativistic Effects 

Beaming. [D’Orazio et al.jpOTS] ) have shown that if PG1302 con¬ 


sists of a massive (M > lO^M©) but unequal-mass (0.03 < 
M 2 /Ml < 0.1) SMBH binary, seen within < 30° of edge-on, 
then the entire 0.14 mag variability of PG 1302 can be explained 
by relativistic Doppler boost. In the hydrodynamical explanation 
proposed here, where the 5.2-yr modulation arises from variations 
in the accretion rate, relativistic boost would also inevitably im¬ 
print additional sinusoidal modulations at the true (shorter) binary 
period. The effect would be enhanced, because the secondary’s ve¬ 
locity is higher by a factor of ^ potentially causing a detectable 
second peak in the periodogram at 5.2x~^ yr. Requiring consis¬ 
tency with § |3.1| we use the maximum binary mass and minimum 
mass ratio {q = 0.03) to put an upper limit on the secondary’s l.o.s. 
velocity uios- The relativistic beaming factor is [r(l — v\os /, 
where V is the Lorentz factor. [P^Orazio et ^ ( |20I5| | have estimated 
the spectral index a = 1.1 from an average over the continuum in 
the V band. We find that the corresponding maximum velocity im¬ 
prints a 0.07 mag amplitude modulation on PG 1302’s light curve. 
PG 1302’s periodogram does not show a significant secondary peak 
with sub-5.2 yr periods, but noise modelling suggests that such sec¬ 
ond peaks would be detectable only at amplitudes of > O.OTmag, 
~ half of the 5.2-yr modulation ( jCharisi et al.|20I5] see the next 
section). 

Iron K a lines. Because the binary separation can be reduced be¬ 
low < lOOi^s, FeK a lines generated at such small separations can 
have characteristic binary-related features, such as ‘missing wings’ 
(due to the central cavity), or ‘see-saw oscillations’ of the red and 
blue wings (due to Doppler-shifting of the emission from mini¬ 
discs; [McKernan et al.||20I3jl. These may be dete ctable with the 
upcoming Astro-H mission { [Takahashi et al.|20I4} . 

3.3 Orbital time-scale Variability 

The binary-i-CBD model discussed above generically predicts mul¬ 
tiple periodic variations. If the observed period of PG 1302 is the 
true binary period, then its periodogram could contain lower fre¬ 
quency, higher-amplitude, and also higher frequency, lower ampli¬ 
tude peaks. These could be revealed in future data, combined with 
more sophisticated search algorithms for periodicity (e.g. jVander-j 
jPlas & Ivezic|2015| and references therein). It will be helpful in such 
a search that two of the periodicities occur at tbin and 0.5tbin, 
i.e. with a characteristic 1:2 ratio in frequency. These are expected 
to be the only two peaks present for 0.05 The variabil¬ 

ity at tbin can disappear entirely, but this happens only in the limit 
ofq^l (presumably rarely realized in nature). Thus the detection 
of a secondary peaks, and the characterization of the full variabil¬ 
ity structure, can help confirm the binary nature of PG 1302, and 
constrain its parameters. jCharisi et al.| ( [20T5] ) searched PG 1302’s 
available photometric data for the existence of additional peaks at 
frequencies above or below the strongest and unambiguous 5.2 yr 
period. No significant peaks were detected, and an upper limit of 
Sm > 0.07 — 0.14 mag (depending on frequency) was derived for 
the amplitude of additional modulations. 

4 CONCLUSIONS 

For binaries with mass ratio in the range 0.3<q< 0.8, hydro- 
dynamical simulations of CBDs predict dominant luminosity vari¬ 
ations at 3—8 times the binary orbital period, due to a dense lump in 
the CBD (Fig.[^. If the periodic variability observed in quasar PG 
1302 is identified with this lump period, rather than the orbital pe¬ 
riod of a putative SMBHB, a two to four times smaller binary sepa¬ 
ration is inferred. This would place the PG 1302 binary securely in 
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the GW-driven regime, making it the first EM detection of such a 
system, and proving that gas can follow the binary past decoupling. 
This is encouraging for the possibility of locating EM counterparts 
of GW sources. Because binaries spend less time at smaller separa¬ 
tions, a shorter tbin is in better agreement with the small number of 
SMBHB candidates reported by G15. The higher orbital velocity 
of the binary increases the effects of relativistic beaming, causing 
optical variability at the orbital period, and also on inferred broad 
line widths. 

The binary-i-CBD model can be tested as it predicts variabil¬ 
ity at multiple, well-defined frequencies which depend on binary 
mass ratio and disc parameters. Since a recent search ( [Charisi et al.| 
|2015| ) did not reveal secondary variability in the optical light curve 
of PG 1302, follow up observations are required. Einally, associat¬ 
ing the BLR with the inner annuli of a lumpy CBD, we find that the 
EWHM of the lines can vary at the period of the continuum vari¬ 
ability by ±14 per cent; we also predict a much smaller shift of the 
broad line centroids. These predictions are consistent with existing 
observations of the width and offset of the broad line. Eollow- 
up spectra, sampling PG 1302’s apparent 5.2 yr period, could test 
this interpretation of the BLR and aid in identifying the nature of 
PG 1302’s variability. 
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